Search for Neutrinoless Double-Beta Decay in 136 Xe with EXO-200 
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We report on a search for neutrinoless double-beta decay of 136 Xe with EXO-200. No signal is 
observed for an exposure of 32.5kg-yr, with a background of ~ 1.5 x 10~ 3 kg _1 yr _1 keV -1 in the 
±lcr region of interest. This sets a lower limit on the half-life of the neutrinoless double-beta decay 
Xe) > 1.6 x 10 25 yr (90% C.L.), corresponding to effective Majorana masses of less than 



r>0 r/gfffl38 
'1/2 



140-380 meV, depending on the matrix element calculation. 
PACS numbers: 23.40.-s, 21.10.Tg, 14.60.Pq, 27.60.+j 
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Double-beta decay is a rare nuclear decay observable 
in those even-even nuclides where /3 decay is either en- 
ergetically forbidden or highly spin suppressed. The 
hypothetical neutrinoless (OV/3/3) decay mode can only 
occur for massive Majorana neutrinos This decay 
may be mediated by the exchange of a Majorana neu- 
trino or by other new particles. The recent discovery 
of neutrino mass in oscillation experiments makes 
the search for the Majorana nature of neutrinos par- 
ticularly relevant and timely. The OV/3/3 decay rate is 
related to the square of an effective Majorana neutrino 
mass {m)/3p by the product of phase space and a nu- 
clear matrix element squared. Kinematic measurements 
restrict the neutrino mass scale to be below 0(1 eV) 0, 
leading to OV/3/3 half lives beyond 10 24 yr. A search 
for OV/3/3 decay in 76 Ge has claimed a positive observa- 
tion [J with T 1 ° / f^( 76 Ge)= (2.23 4 :°; 4 t) x !0 25 yr, implying 
( m )/3/3 = 0.32 ± 0.03 cV for the nuclear matrix element 



given in 

The exceedingly long half lives of interest for 0vf3(3 
decay require large detectors using isotopically enriched 
sources, radio-clean construction techniques and the abil- 
ity to actively reject remaining backgrounds. The Oi//3/3 
decay results in a discrete electron sum energy distribu- 
tion centered at the Q- value (Qpp). The allowed, yet also 
rare, two neutrino double beta (2vfif3) decay is character- 
ized by a continuous sum energy spectrum ending at i 
The 2v(3f3 decay has been observed in many isotopes _ 
and, recently, in 136 Xe with the EXO-200 detector 
later confirmed in Q. The decay electron sum energy 
allows discrimination between the two modes. 

EXO-200, described in detail in Q, uses xenon both 
as source and detector for the two electrons emitted in 
its decay. The detector is a cylindrical homogeneous 
time projection chamber (TPC) |9|. It is filled with liq- 
uefied xenon ( enr LXe) enriched to (80.6±0.1)% in the 
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isotope 136 Xe. The remaining 19.4% is 3 Xe, with other 
isotopes present only at low concentration. EXO-200 is 
designed to minimize radioactive backgrounds, maximize 
the enr LXe fiducial volume, and provide good energy res- 
olution at the 136 Xe Q- value of 2457.83 ± 0.37 keV 0. 
Energy depositions in the TPC produce both ionization 
and scintillation signals. The TPC configuration allows 
for three-dimensional topological and temporal recon- 
struction of individual energy depositions. This ability 
is essential for discriminating /3(3 decays from residual 
backgrounds dominated by 7s. 

The cylindrical TPC is divided into two symmetric 
volumes separated by a cathode grid. Each end of the 
TPC is instrumented with 38 charge induction (V) and 
38 charge collection (U) wire triplets. The U and V wire 
grids, crossing at 60°, provide stereoscopic information 
for charge depositions. At each end of the TPC there 
are ^250 Large Area Avalanche Photodiodes (LAAPDs) 
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that record the 178 nm scintillation light. A drift 
field of 376 V/cm is applied in the TPC volume. All sig- 
nals are digitized at lMS/s. All detector components 
were carefully selected to minimize internal radioactiv- 
ity [Hj]. The TPC is mounted in the center of a low- 
background cryostat. At least 50 cm (25 cm) of HFE-7000 
fluid [13[ (lead) shield the TPC from external radioactiv- 
ity. Calibration sources can be inserted at various posi- 
tions immediately outside of the TPC. The clean room 
module housing the TPC is surrounded on four sides by 
50 mm thick plastic-scintillator cosmic-ray veto panels, 
which are (95.5 ±0.6)% efficient. EXO-200 is located at 
a depth of 15851" m.w.e. at the WIPP in New Mexico, 
USA. 

EXO-200 started taking low background data in late 
May 2011. The data presented here were collected from 
September 22, 2011 to April 15, 2012, for a total of 
2,896.6 hours live time under low background conditions. 
During the same period, 376.8 hours of calibration data 
were collected with three 7 sources at three positions. 
Most of these data was taken with a 228 Th source for the 
primary purpose of measuring the free electron lifetime 
(r e ) in the TPC. This is achieved by recording the charge 
collected for 2615 keV full absorption events occurring at 
different z-positions. The resulting r e was ~200 ^s for the 
summer 2011 but climbed to ~3ms upon increasing 
the Xe recirculation flow for the current dataset. 

Offline event reconstruction proceeds in three stages: 
signal finding, parameter estimation, and clustering. 
Charge signals on the U wires and scintillation signals 
on the two LAAPD planes are found using a matched 
filter technique. The filter yields time estimates for both 
ionization and scintillation channels. The time informa- 
tion from the U wires is used to search for induction 
signals in V wire waveforms. A signal un-shaping algo- 
rithm is applied to the U wire waveforms, optimizing the 
discrimination between single-site (SS) and multiple-site 
(MS) energy depositions. Candidate U wire signals are 
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FIG. 1: Correlation between ionization and scintillation for 
SS events from a 228 Th source. The energy resolution is con- 
siderably improved by forming the linear combination of both 
measurements. Events in the top-left quadrant are due to in- 
complete charge collection and are rejected by the cut (dashed 
line), removing only 0.5% of the total. 



fit to template waveforms modeling the measured trans- 
fer functions so that the signal amplitudes can be ex- 
tracted for energy estimation. Amplitudes are corrected, 
channel-by-channel, for electronic gains determined from 
radioactive source calibration. U wire and V wire sig- 
nals are then combined into charge clusters using timing 
information from the fits, and associated with the near- 
est (earlier in time) summed scintillation signal. Each 
cluster energy is corrected for position-dependent charge 
losses due to finite enr LXe purity and for the shielding 
grid inefficiency of the V wire plane. This procedure, in- 
troducing the only time-dependent correction, yields a re- 
constructed ionization energy and three dimensional po- 
sition information for each charge cluster. An efficiency 
loss is incurred by events for which 3D reconstruction is 
not possible as these are rejected from this analysis. 

The LAAPD sum signal is corrected for variations 
of gain and spatial light collection efficiency within the 
TPC. A trilinear interpolation of the 2615 keV 7 SS full 
absorption peak scintillation signal, recorded for 1352 lo- 
cations throughout the 6nr LXe during 228 Th calibration 
source deployments, is used to construct the correction 
function. 

Figure Q] shows the energy of events as measured by 
the ionization and scintillation channels while the 228 Th 



source was deployed. As first discussed in [14j and evi- 
dent from the tilt of the 2615 keV full absorption ellipse 
in the figure, the magnitude of the two signals is anti- 
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FIG. 2: Top: systematic uncertainty bands on the energy 
calibration residuals, using the full energy reconstruction de- 
scribed for the three 7 sources. For both SS (solid) and MS 
(dashed) the position of the four 7 lines is consistent with the 
calibration model within < 0.1%. Bottom: energy resolution 
for various sources along with a fit to the empirical model (see 
text). 



correlated. The 2D SS and MS energy spectra are in- 
dependently rotated and projected onto a new (ID) en- 
ergy variable in such a way as to minimize the width of 
the 2615 keV 7 line. Energy spectra from 137 Cs, 60 Co, 
and 228 Th sources are produced using this method and 
then, the positions of the full absorption peaks at 662, 
1173, 1333 and 2615 keV are fit. The energy calibration 
function that converts the rotated energy estimator into 
keV has a small quadratic term. The residuals, defined 
as (Eftt — E tlne ) / Et rn e, and the uncertainty bands are 
shown in Figure [2] (top) . The energy resolution, parame- 
terized as cr 2 = aa 2 +bE+cE 2 , is plotted versus energy in 
Figure [2] (bottom). Here a e is the electronic noise contri- 
bution, bE represents statistical fluctuations in the ion- 
ization and scintillation, and cE 2 is a position- and time- 
dependent broadening. It is found that a/E = 1.67% 
(1.84%) for SS (MS) spectra at Qpp, dominated by the 
noise and broadening terms. 

The ability of the TPC to identify SS and MS inter- 
actions is used to separate (3 and /3/3 decays in the bulk 
xenon from multiple site 7 interactions. The clustering, 
currently applied in 2D, has a separation resolution of 
18mm in the U-dimcnsion and 6mm in z (drift time). 
The event sharing between SS and MS energy spectra 
is demonstrated for a 228 Th source in Figure |3l The 
SS and MS spectra are compared to probability density 
functions (PDFs) generated by GEANT4 [H] simulations 
(MC). Simulated charge depositions are transported to 
the wires where signals are generated using a model of 
the readout electronics. Noise is added and signals are re- 
constructed in the same manner as for data. After events 
are designated as either SS or MS, the total collected 
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FIG. 3: MS (top) and SS (bottom) energy spectra from a 
228 Th calibration run. The PDF generated from simulation 
(line) is fit to the data floating the normalization of the MS 
and SS spectra independently to illustrate the agreement with 
the rendered spectral shape. The limited ability of the model 
to reproduce the absolute source rate is quantified in the text 
and is not illustrated by the figure. 



energy is convolved with the energy resolution function 
shown in Figure [2] This procedure reproduces the shape 
of the 60 Co and 228 Th source spectra (see Figure [3] for 
228 Th). The MC reproduces the fraction of SS events, 
defined as N SS /(N SS + N MS ), to ±8.5%. In addition, for 
228 Th (Qu/3P) events the simulation predicts a 70% (71%) 
efficiency for the requirement that events are fully re- 
constructed in 3D. The simulation predicts the absolute, 
NIST-traceable, activity of all sources to within ±9.4%. 
These variations are used as measures for the systematic 
uncertainties in the final fit. The 71% MC estimate for 
the 0vf3f3 efficiency is further verified over a broad en- 
ergy range by comparing the 2i//3/3 MC efficiency with 
low-background data. Although the 2i//3/3 spectrum has 
vanishing statistical power at Qpp, the efficiency is found 
to be a smooth function of the energy and agrees with 
the simulated efficiency to within the ±9.4% overall scale 
uncertainty mentioned above. 

The fiducial volume used in this analysis contains 
79.4kg of 136 Xe (3.52xl0 26 atoms), corresponding to 
98.5 kg of active enr LXe. The trigger is fully efficient 
above 700 keV. The cut represented by the dashed diago- 
nal line in Figure [1] eliminates a population of events due 
to interactions in the enr LXe region for which the charge 
collection efficiency is low, leading to an anomalous light- 
to-chargc ratio. This cut also eliminates a decays from 
the low background data, but causes only a negligible 
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FIG. 4: MS (top) and SS (bottom) energy spectra. The 
best fit line (solid blue) is shown. The background com- 
ponents are 2i//3/3 
(dotted dark blue) 



grey region), K (dotted orange), Co 
!22 Rn in the cryostat-lead air-gap (long- 



dashed green), i3S U in the TPC vessel (dotted black), 2d2 Th 
in the TPC vessel (dotted magenta), 214 Bi on the cathode 
(long-dashed cyan), 222 Rn outside of the field cage (dotted 
dark cyan), 222 Rn in active xenon (long-dashed brown), 135 Xe 
(long-dashed blue) and 54 Mn (dotted brown). The last bin on 
the right includes overflows (none in the SS spectrum). 
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FIG. 5: Energy spectra in the 136 Xe Qpp region for MS (top) 



and SS (bottom) events. The I (2)cr regions around LJpp are 
shown by solid (dashed) vertical lines. The Ov/3/3 PDF from 
the fit is not visible. The fit results have the same meaning 
as in Figure |4] 



loss of efficiency for 7- and /3-like events. Cosmic-ray in- 
duced backgrounds are removed using three time-based 
cuts. Events preceded by a veto hit within 25 ms are re- 
moved (0.58% dead time). Events occurring within 60s 
after a muon track in the TPC are also eliminated (5.0% 
dead time). Finally, any two events that occur within 1 s 
of each other are removed (3.3% dead time). The combi- 
nation of all three cuts incurs a total dead time of 8.6%. 
The last cut, combined with the requirement that only 
one scintillation event per frame is observed, removes /3- 
a decay coincidences due to the time correlated decay 
of the 222 Rn daughters 214 Bi and 214 Po. Alpha spectro- 
scopic analysis finds 360 ±65 /iBq of 222 Rn in the enr LXe, 
that is constant in time. 

The SS and MS low background spectra are shown in 
Figure [4] Primarily due to bremsstrahlung, a fraction 
of /3/3 events are MS. The MC simulation predicts that 
82.5% of Of/3/3 events are SS. Using a maximum like- 
lihood estimator, the SS and MS spectra are simultane- 
ously fit with PDFs of the 2f /3/3 and Ov/3/3 of 136 Xe along 
with PDFs of various backgrounds. Background models 
were developed for various components of the detector. 
Results of the material screen campaign, conducted dur- 
ing construction, provide the normalization for the mod- 
els. The contributions of the various background com- 
ponents to the Of/3/3 and 2f/3/3 signal regions were esti- 
mated using a previous generation of the detector simula- 
tion For the reported exposure, components found to 
contribute < 0.2 counts (Of/3/3) and < 50 counts (2f/3/3), 
respectively, were not included in the fit. For the current 
exposure, the background model treats the activity of the 
222 Rn in the air-gap between the cryostat and the lead 
shielding as a surrogate for all 238 U- like activities exter- 
nal to the cryostat, because of their degenerate spectral 
shapes and/or small contributions. A possible energy off- 
set and the resolution of the 7-like spectra are parameters 
in the fit and are constrained by the results of the source 
calibrations. The fraction of events that are classified 
as SS for each of the 7-like PDFs is constrained within 
±8.5% of the value predicted by MC. This uncertainty 
is set by the largest such deviation measured with the 
source calibration spectra. The SS fractions for /3- and 
/3/3-like events are also constrained in the fit to within 
±8.5% of the MC predicted value. As a cross-check, the 
constraint on the 2f/3/3 SS fraction is released in a sep- 
arate fit of the low background data. The SS fraction is 
found to agree within 5.8% of the value predicted by the 
MC simulation. 

The /3/3 energy scale is a free parameter in the fit, so 
that it is constrained by the 2f /3/3 spectrum. The fit re- 
ports a scale factor of 0.995 ±0.004. The uncertainty is 
inflated to ± 0.006 as a result of an independent study of 
the possible energy scale differences between 7- and /3/3- 
like energy deposits. The 2f/3/3 PDF is produced using 
the Fermi function calculation given in [16| . Tests using 
a slightly different spectral form (l7j were performed and 
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constraint on this cosmogcnic background, the smaller 
limit is reported. 

The result from the likelihood fit is shown in Figure |6l 
along with the recent constraint for 136 Xe Q and the 
best limit [lj| and claimed detection [3] for 76 Ge. The 
present result contradicts [3] at 68% C.L. (90% C.L.) for 
the nominal values of all (most) matrix element calcula- 
tions considered 

muni! 

and provides upper bounds to 
Majorana neutrino masses between 140 and 380 meV at 
90% C.L.. 
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Relation between the T^ p in 76 Ge and 136 Xe for 
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different matrix element calculations (GCM [H, NSM [HJ 
IBM-2 [12|, RQRPA-1 [H| and QRPA-2 H). For each matrix 
element (m)pp is also shown (eV). The claim Q is represented 
by the grey band, along with the best limit for 76 Ge fTgf| . The 
result reported here is shown along with that from 



found to contribute <0.001% to the predicted location 
of Qpp- Finally, the stability of this energy scale cor- 
rection was verified by repeating the fit with increasing 
thresholds up to 1200 kcV. 

For the best-fit energy scale and resolution the ±ler 
and ±2ct regions around Qpp are shown in Figure[5] The 
number of events observed in the SS spectrum are 1 and 
5, respectively, with the 5 events in the ±2<r region ac- 
cumulating at both edges of the interval. Therefore, no 
evidence for Oi//3/3 decay is found in the present data set. 
The lower limit on T®J^ is obtained by the profile likeli- 
hood fit to the entire SS and MS spectra. Systematic un- 
certainties arc incorporated as constrained nuisance pa- 
rameters. The fit yields an estimate of 4.1 ± 0.3 back- 
ground counts in the ±ler region, giving an expected 
background rate of (1.5 ± 0.1) x 10~ 3 kg _1 yr _1 keV~ 1 . 
It also reports 0^/3/3 decay limits of < 2.8 counts at 
90% C.L. (< 1.1 at 68% C.L.). This corresponds to a 



T°f p > 1.6 x 10 25 yr at 90% C.L. (T^ p > 4.6 x 10 25 yr 
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at 68% C.L.). Toy MC studies confirm the coverage of 
this method as suggested by [HI- The same fit also re- 
ports T^ 13 = (2.23 ± 0.017 stat. ± 0.22 sys.) x 10 21 yr, 
in agreement with @ and 0. The levels of contamina- 
tion from 7-emitting nuclides are found to be consistent 
with material screening estimates 12J. The addition to 
the fit of a PDF for 137 Xe produces a 13% higher limit 
on Ty^P at 90% C.L.. In the absence of an independent 
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